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I discuss and extend the recent proposal of Leclair and Mussardo for finite temperature 
correlation functions in integrable QFTs. I give further justification for its validity in the 
case of one point functions of conserved quantities. I also argue that the proposal is not 
correct for two (and higher) point functions, and give some counterexamples to justify that 
claim. 
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In a recent paper [PJ, A. Leclair and G. Mussardo have come up with a very interesting 
suggestion to compute correlation functions in integrable quantum field theories at finite 
temperature. In a nutshell, what they propose is to identify the usual particle and hole 
excitations over the "thermal ground state" || , with their non trivial dressed energies and 
momenta (as obtained by solving the thermodynamic Bethe ansatz equations), but still 
take, for the physical operators, the undressed form factors. 

The proposal in 0] is quite exciting, as finite temperature correlation functions are 
quantities of the highest interest (in particular for comparison with experiments). However, 
it is fair to say that no real demonstration is provided in except for some one point 
functions, and that the formula in JIJ are merely a tentalizing guess. 

I will argue in this note that, while the formula proposed in |1| for one point functions 
are, at least for some operators (the densities of conserved quantities), correct, the ones 
for two and higher point functions generally do not hold. The reason for this is quite 
simple and physical: the formula proposed in |IJ does not take into account the effects 
of the dressing in the form-factors. This is fine for one point functions: since the same 
multiparticle states are on the left and on the right side of the correlator, there is no motion 
of the thermal ground state due to interactions, and the formula for bare form factors 
can still be used. In contrast, for two point functions, one has to consider intermediate 
states where particle or holes have been created; this leads to a displacement of the thermal 
ground state, and a dressing of the form- factors, that has to be taken into account - there 
is just no reason for the bare form factors to still be relevant there. 

One of the difficulties in assessing the validity of [0 is that very few finite temperature 
correlators are known exactly in integrable theories, and that the expressions proposed in 
0] , even if conceptually quite simple, are nevertheless very hard to compute explicitely 0. 
However, as I will argue below, the argument presented in [|l| is quite general, and holds 
just as well for correlators evaluated in other thermodynamic ensembles. The case of a 
chemical potential at T = turns out to be particularly simple, and will allow me to put 
the two point functions of to a serious test, that they unfortunately fail. 

Let me now proceed and discuss one point functions first. I will only consider operators 
0(x) (I refer to them, a bit incorrectly, as conserved quantities) for which the quantity 
O = J 0(x)dx acts diagonally on multiparticle states, with one particle eigenvalues o(6). 

1 A possibility might be to investigate the low temperature limit, and compare it with results 
obtained in Q , but I won't do this here. 
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Examples of this include the energy and momentum operators, and after a very slight 
modification, the trace of the stress energy tensor considered in By translational 
invariance, it thus follows that 



(e n ,...,e 1 \o{x)\e 1 ,...,e n ) = iA 

{9 n , . . . , 9i\9i, . . . , 9 n ) L 



i=l 

L the system length. 

I want now to compute the average of O at finite temperature, and in the presence 
of possible other thermodynamic couplings. I will restrict for notational simplicity to a 
theory with a single particle, and consider, in addition to the temperature, the presence 
of a chemical potential. Generalizations are quite straightforward. Due to the non trivial 
S matrix of the theory, the density of allowed states P (P = p + p h ) obeys (I set ft = 1 
here) : 

2nP(9) = Lm cosh # + 2tt$ * (1) 

where I defined / * g(9) = J f(9 - 9')g{9')^. In ©, the kernel $ = In 5, where S 

is the scattering matrix. 

The average of O is computed using the thermodynamic Bethe ansatz method: con- 



figurations are weighed with a Boltzmann weight exp 



E-TS-nN 



where T is the tem- 



T 

perature, S the entropy, p the chemical potential, and iV the number of particles, and 
saddle point equations are written, leading to the TBA equations. Introducing the pseudo 

h , 

energy e defined by — = e e , the TBA equation at temperature T ( T = j^) read 

e(0) = mR cosh#- $*ln(l + e" e+Mi? ) , (2) 

It can be shown with the usual argument [|J that 2irP = L §^\ flR f ixed - From this the 
average of O follows: 

(0(x)) =\ [ o(9)p(9)d9 



L 



d9_ 
2^' 



jiR fixed 

with /_ = 1+e e„ M R ■ I can now solve the TBA equation iteratively to obtain the expansion 



(0(x))= [ o{9) f-(9)m cosh 9—+ [ o(9)f_(9)— ! $(9 - 9') f_(0') m cosh 9' — 
J 27T J v 2tv J 2tv 

o(9)f-(9)^ J $(6>-0')/-(0')^ / $(9'-9")f-(9")mcosh9"^ + ... 



Notice that this expression in terms of o, $ and the filling fraction /_ is completely general. 
It would hold for a system with more thermodynamic couplings, provided these couplings 
involve conserved quantities. The expression (|J) also generalizes easily to cases with several 
types of particles. 

I will now OD) compare with the formula proposed in [0: 

<°( X )> = E^ / ^■■■^f[f-(0i)(0n...9i\o(x)e 1 ...e n ) conn , (5) 

n=0 ' ^ i=l 

and argue that it is the same. Of course, in |IJ, formula (^j) is written only in the case of 
finite temperature. However, all arguments presented there rely on cancellations between 
numerators and denominators, together with hints from the free case and the one point 
function of the stress energy tensor, and these are all features which generalize to more 
complicated thermodynamic averages, like the ones involving a chemical potential. In (||), 
the connected form factor is obtained in principle by using the crossing formula for non 
coinciding arguments, and then getting rid of all the diverging or ill defined terms as the 
arguments are sent to one another. 

The key point here is that the values of the connected form-factors can be understood 
generally as follows (this argument is quite similar to Balog's construction ||). To start, 
I need to make some remarks on the normalization of states. Consider the scalar product 
(9\9), which, formally, is equal to 2nS(0), since the normalization used in the construction 
of the asymptotic states is (9\9') = 2tt5(9 — 9'). Of course the symbol 5(0) does not make 
much sense, and has to be regularized properly. The way to do this is to consider the 
completude relation, 1 = J ^-\9)(9\. This relation hides the cancellation of two terms, the 
meaning of which is easier to see by introducing a finite length L in the system. First, the 
state \9) is actually not normalized, because of the 5(0) term mentioned just before; the 
integral should therefore involve instead the ket )— — and similarly for the associated bra. 
Second, allowed rapidities for a particle are not uniformly distributed; rather, they have a 
density given by 2txP = mL cosh 9, and the completude relation should involve therefore 
an integral J P(9)d9. For the completude relation to be equivalent to the one we used so 
far, which preserves the scalar product of states, we thus need to have || 

(9\9) =mL cosh 9. (6) 
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The argument can be generalized to the case with several rapidities. For two particles, 
the Bet he equations read: 



2-7mi =mL sinh0i H — \nS(9i — 62) 

1 

2tvu2 =mL sinh 6*2 H — In S(02 — 61 ) ■ 

From this, we deduce that 

mL cosh 6>i + $(0i 2 ) - $(0i2) 

(0 2 ,0i|0i,0 2 >=Det 

-$(0 i2 ) mL cosh 2 + $(0i 2 ) 
=m 2 L 2 cosh0i cosh 02 + mL$(0i2)(cosh0i + cosh 02). 



(7) 



(8) 



The formula generalizes to an arbitrary number of rapidities. As could have been expected, 
it coincides with the well known formula for the norm of Bethe states derived by Gaudin 
and by Korepin in the context of the XXZ and other models (see |7j). 

I can now get back to @ and (]5|). It is clear that the first term of the expansions 
match, due to (§) 

(6\0(x)\6) = (0|O(x)|0) corm = o(0)mcosh0. (9) 



For the two particle form factor, let me define: 

(02,0l|O(x)|0l,0 2 > C onn =<0 2 , 01 \0 (x) 1 01, 2 ) - (0 2 \0(x) 1 2 ) conn (6 X \ 6 X ) 

-(0i|O(x)|0i) conn (0 2 |0 2 ). 



(10) 



Here, the scalar products (eg (0i|0i)) have to be computed in the presence of the other 
particle (02), something I left implicit in the notations for simplicity. From the general 
arguments explained previously, (0i|0i) = mL cosh 0i + </>(0i2)- It thus follows that 

(02, 01 |0(X) |0l, 2 ) CO nn =m$(0i 2 ) [cOSh 01O(0 2 )+ COSh 2 O(01)], (11) 

again ensuring a matching between (|j) and ((^). A general formula for connected form- 
factors follows easily by extending (|T0| ) from n = 2 to arbitrary values of n: simply subtract 
from the initial form factor all the possible contractions, being careful to evaluate scalar 
products of states in the presence of all other particles - they are all expressed as various 
minors of the same initial determinant. The net result is simply: 

(0 n , . . . , 0i|0(o;)|0i, . . ..On) conn = m$(0i 2 )$(0 2 3) • • • $(0n-i,n) cosh0io(0 n )+ permutations, 

(12) 
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and clearly ensures the coincidence of (Q) and (|5|). 

I am of course not pretending that this argument is a proof of ([5]), but I believe it 
could be sufficiently strenghtened, maybe along the lines of ||. As it stands, it certainly 
gives further support to ([5]), and to my claim that the guess of should be considered 
for more general thermodynamic ensembles. 

I shall now argue that, unfortunately, the formula proposed in |l| for the two particle 
correlator is probably not correct. Rather than give general arguments, I will present 
a simple counterexample, in the case where T = but there is a non trivial chemical 
potential. 

Consider indeed the free boson with hamiltonian 



2 



8ngU 2 + — - (8 X 
8ng 



dx. (13) 



I will add to this hamiltonian a term J d x (pdx. Properties in the presence of this coupling 
are straightforward to evaluate by a simple shift of the bosonic field. Introducing chiral 
components 99, (p, one finds (d x 4> = d z (p + dz'fi) 



(d x <p) =Vg 
(d Zl <pd Z2 <p) =(Vg) 



2g (14) 



where I have set z = x + iy, y the imaginary time (called t in [|l]). 

On the other hand, we can consider the free boson as the UV limit of the sine- Gordon 
theory, where the hamiltonian ([13|) is supplemented by an interaction term A j cos<pdx. 
Provided we consider physics at a scale much smaller than the correlation length induced 
by this perturbation, oc A _2+2a , we will observe results similar to the free boson case. 
Equivalently, we expect to be able describe the properties of the free boson theory using 
a massless scattering description [Q, with purely right and left moving particles obtained 
by taking the large rapidity limit of the usual solitons, antisolitons and breathers of the 
theory. The parametrization of the energy I will use in that limit is e = ±p = e 9 , where 
I have set an arbitrary mass scale equal to unity. The correlators (|14|) should therefore 
be obtainable using the formulas proposed in 0. I will show that, actually, only the one 
point function is obtained. 

For simplicity, I restrict to the case T = 0. The field V then leads to the creation 
of a pair of Fermi seas of left and right moving massless solitons: in the following, I will 
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concentrate on the right moving sector only. The average of d z ip is directly related with 
the number of solitons in the ground state 



-L(J dz pdx) = (N) 



This quantity fits exactly in the context described previously. The two particle kernel 
is given by the soliton soliton S matrix, and, since we are at temperature T = 0, the 
function /_ becomes very simple: /_ = 1 for particles in the sea, /_ = otherwise. If 
we consider the operator d z (p, it couples only to right moving particles, for which the sea 
is G (— oo, A), A a Fermi rapidity. The average of d z <p follows from the appropriate 
generalization of ([|): 

Expressions for the form- factors of the operator d z tp are well known: we can, in 
principle, directly compute the connected ones, and show agreement with the general 
formula (|12]). For instance, the two particle form- factor gives 

(9\d z <p\6) = 2ne\ (16) 

in agreement with (the massless limit of) @ for o(0) = 27r, the value of the integral of 
d z ip for a soliton. I have checked similarly the formula for the connected four particle 
form- factor (|TT|) for simple values of the coupling g. 

I now turn to the two point function of the operator d z ip, and show that the expression 
proposed for it in [[TJ this time does not work. Recall the proposal of M: 

(O(x,y)O( ,0)) R =((O) R f + ±± £ .^fUft) 

AT=1 ' (Ji=±l J j=l (17) 

exp (i(jj{xkj + iyej)^j |(O|0(O)|0i, . . ■ ,0 N ) ai: „, )aN \ 2 . 

Here, /+ = 1+ ^,_ £ , e = Re is the dressed excitation energy, while p is the dressed excitation 
momentum. Finally, 

(010(0)10!, . . . , 0^,...,^ = (O|O(O)|0! - m(a 1 - l)/2, ...,6 N - in(a N - l)/2> 

The physical meaning of (|I7]) is more transparent than the formula. Correlations at temper- 
ature T should be determined by processes involving excitations over the thermal ground 
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state, with dressed energy and momentum 0. These excitations can be of two types, 
particles or holes (the variable a in (|PTD ). A sum as fll7|) is thus expected; the key proposal 
of is that the form factors for the physical excitations are the same as the bare ones. 

In the case we consider, the situation simplifies considerably. First, let me recall the 
structure of the excitations B. It is conveniently represented in the following table: 



holes : e\(6), 9 < A 

solitons : e+(0), 9 > A 

antisolitons : £-(0) > (1 — g)V, 9 arbitrary 

breathers : e n (9) > ngV, 9 arbitrary 



(18) 



The excitation energies e have non trivial expansions given in 0. 

The physical processes involved in the two point function of d z (p are of three basic 
types: creation of a pair particle hole (ie take a soliton in the Fermi sea and move it 
outside); creation of a pair soliton (above the Fermi sea) antisoliton, and creation of a 
breather. The latter two processes have thresholds. Since the quantities e and p have 
no singularity, and the bare form-factors do not know anything about these thresholds, 
it is immediately quite obvious that the proposal of |1|] cannot be true: it would lead to 
singularities in the Fourier spectrum of the correlator, in sharp contrast with the result 
expected from ([14]). 

Nevertheless, it is probably reasonable to make this counter example more explicit. 
First, let me show that the formula is in fact right for the free case g = \. In that case, 
things simplify for several reasons: there are no breathers, the excitation energies have 
simple expressions e\ = ¥ — e e , e+ = e 9 — ¥, e_ = e e + \ , and only the two particle form 
factor of d z (p is non zero, (O\d z (p\0i02) = 2iTce 9l ^ 2 e e2 ^ 2 . The connected term in expression 
([17D reads then 



+ 



d9 



r°° 

' J A 


d9 2 
2tv 


2ne 




f°°d9 2 


J 


A 27T 



A/V 2 / 2 



h/2 e 9 2 /2 



1 



— [l - exp(e A z)] j exp(e 



A, 
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the expected result ([14]) for g = 

In the general case, I will consider the Fourier transform of the two point function, 
S(u) = J ^e~ luJX (d z f(x,y)d z if(0)) y ^o = 2gu {u > 0). I will also restrict to the case of 



low frequencies (that is, ui smaller than any of the thresholds), where only the particle hole 
processes contribute; if formula (|17|) fails in that case, it will be enough to show it is not 
correct. 

Below the lowest thresholds, S(u), according to formula (0), is given by a sum of 
terms S n corresponding to processes with n particles and n holes. When V = and we 
know formula (jTj]) is correct, each of the S^s is linear in u, S n = c n u. In that case, there 
are no thresholds, so the S n have to be supplemented by the terms corresponding to the 
other processes, each of which is also proportional to 10. The infinite sum of all these ui 
terms reproduces the desired behaviour S(u>) = 2gu; as checked in |TTJ the convergence of 
this sum is in fact very quick for g not too close to 1. When V ^ 0, formula ( |TT| ) predicts 
correctly (this simply follows from dimensional analysis) that all the <S n 's now have the 
form S n = uf n (y), with c n = f n (oo). Let us now consider the limit of small frequencies 
at finite V: in that case, we are instead exploring the behaviour of the functions f n in 
the limit of very small argument. Since the contribution at frequency uj is determined by 
rapidities 6i > A of holes and 6[ < A of particles such that 

n n 
i=l i=l 

clearly, for small u, the particles and holes have to be closer and closer to the Fermi 
rapidity. The two particle contribution for instance, which reads in general 

f A f 1 r^\(e[\dMei)\ 2 t["-z h + (0i)-t + (e'i)} 

J-oo 2n J A 2n 

becomes in the limit of small u, \ (Ad z (p\A)\ 2 ^ + ^ k h py| (dots denote derivatives with 
respect to the rapidity variable). In fact, for the limit lim^^o ^77^ j this two particle 
contribution is the only one to consider. There are two reasons for that: one is that the 
next contribution is down by an w 2 term due to phase space considerations; the other is 
that in the limit u — >• 0, the rapidities of the n particles and n holes are all compressed 
towards A, and form- factors like eg (A, A\d z (p\A, A) vanish due to the general behaviour 
under rapidity exchanges (S(0) = —1). It follows that, according to [Q, one would find 



The right hand side can be evaluated using the detailed results of || for the filling of the 
ground state etc. One finds (for dimensional reasons, this has to be a pure, V independent 
number) 

RHS = L^(W- S )l\\ (21) 
2*1-9 Vr[l/2(l-g)]J 

Except when g = \, the right hand side is off the exact result (it should be equal to 2g) 
by a finite amount, demonstrating that ([17]) is not correct in that case. 

The reader might worry here that the whole argument assumes somehow absolute 
convergence of the series in (fi7|), and that maybe one would get the correct limit in (|2C| ) 
by first summing over all contributions, then letting uj — > 0. The answer to this is that I 
have considered the limit (|20|) only to make things as clear as possible. We are of course 
interested not only in this limit, but in the behaviour of S{u) on the whole real axis, where 
we have to recover S(u) = 2gu for fllTf ) to be correct. As argued above, only the one 
particle hole contribution has a term linear in uj: other terms start with higher powers in 
u, and it is easy to see that only a finite number of them contribute to a given order in 
u. The series representing S(u) according to (|T?|) is thus of the form ^2dkU k , with all 
the dfc's finite (eqn. (^1]) means that d\ is not the right one; one can also check that d% is 
not the right one - that is, it does not vanish - etc). For my argument to be spoilt by a 
convergence problem, one would need this series to diverge and to somewhat "represent" 
the simple linear term 2gu; besides the fact that I have found no indication of divergence 
numerically (that is, convergence seems as good as in ||1Q|1 ), I do not think such a scenario 
is likely at all. 

The reader might also be surprised by the fact that (|17|) is right for V = but not 
V 7^ 0. The point is once again that ( [P7D leads to a representation of S(u) as a sum of 
terms of the form uf n (y); validity in the case V = is a statement about the sum of 
/ n (oo)'s, while invalidity in the case V ^ is a statement about the general shape of the 
sum of / n 's at finite argument. 

The physical origin of the failure of ( [PTD is easy to trace back to the dressing effects. 
In fact, in 0, another approach to compute correlators at T = with a field coupled to the 
£7(1) current was proposed. In this paper, it is recognized that dressed excitations must 
have dressed form-factors; an expression was proposed in particular for the low energy 
behaviour of the particle hole form factor, which reproduced S(u ) correctly. 

I should stress that I have not found simple counterexamples in the context originally 
considered in pj, of a theory with a temperature and no chemical potential. However, I 
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believe the argument presented in p]] generalizes to any thermodynamic average, so this 
counterexample at T = with a chemical potential is nevertheless a good, alas negative, 
test, of their result for two and higher point functions. Some easy things can be said about 
the T^O case however. One of them has to do with thresholds. Consider for instance the 
massless limit of the sine-Gordon model once again, and its massless scattering description. 
At T 7^ 0, it is one of the striking features of the interactions that the energy necessary to 
add a particle has a gap: for instance, the energy to create a one-breather in the attractive 
regime is Te\ > Tin 3; similarly the energy gained by destroying such a breather is larger 
or equal to this number (this follows simply from the solution of the TBA equations). As a 
result, if we consider again the correlator of d z ip, processes involving the one-breather have 
thresholds (recall that the bare form- factors of d z ip with an even number of one- breathers 
are zero). Since the form factors in the sum ([17]) are the bare ones, they know nothing 
about these thresholds; as for the other pseudo energies, they have no singularity at the 
position of these thresholds. It follows that, according to ([T7[), the Fourier transform of 
the two point function of d z (p would exhibit singularities at finite values of the frequency 
uj. This is of course in contradiction with the simple form of this two point function , 
that follows from conformal invariance, and indicates once again that the proposal in Jl ] 
is generally not correct. 

In conclusion, I believe that the formulas of [|J for the two and higher point functions 
are in general incorrect because they don't take properly into account the dressing of the 
effective vacuum created by finite thermodynamic couplings. As for the one point function 
of conserved quantities, the formula of |I| looks very reasonable when compared with the 
result of the TBA, and presumably could be rigorously proven by a more serious analysis 
of connected form-factors than the one I have presented here. I am not sure about the 
one point functions of non conserved quantities - the examples studied in |l| are quite 
convincing, but further work is probably needed to settle the issue. 
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